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INTRODUCTION 


Because of many unique and attractive properties (e.g., light weight, low cost, 
and compatibility with integrated circuits), microstrip antennas are finding many 
applications. As a result, in recent years a number of theoretical and experimental 
methods to determine their electromagnetic properties have been proposed (ref. 1). 

By treating the rectangular antenna as two parallel slots interconnected by a low- 
impedance transmission line, it is possible to obtain some useful predictions of the 
radiative properties (refs. 2 and 3). However, this simple model is not adequate for 
predicting the variation of impedance with feed location. Furthermore, the 
transmission-line model cannot be applied to microstrip patch antennas other than 
rectangular. To analyze a wider variety of microstrip patch shapes, a cavity model 
has been proposed (ref. 4). Although the cavity model predicts quite accurately the 
radiation pattern, it requires a semi empirical correction for the resonant frequency 
in order to match the calculated results with the experimental data. This discrep- 
ancy may be due to the fringing field at the edges of the patch not being incorpo- 
rated in the formulation and the dielectric substrate with finite thickness being 
replaced by a homogeneous medium of equivalent dielectric constant. There is, there- 
fore, a need for more rigorous analysis techniques for microstrip antennas. 

The rectangular and circular microstrip antennas have been analyzed by modeling 
the patch as a grid of wires and solving the resulting structure numerically 
(ref. 5) . Although the method appears to be more accurate and general, it requires a 
lot of computer time and storage. Furthermore, to optimize the computation, it is 
essential to know a priori the direction of current flow on the patch. 

Using the reaction integral equation in conjunction with the method of moments, 
Newman et al. (ref. 6) propose a method which seems to be quite accurate for predict- 
ing patch current, impedance, and resonant frequency. However, Newman* s method 
requires an unusually precise computation of the elements in the impedance matrix 
while still modeling the dielectric substrate as an equivalent semi -inf inite homoge- 
neous medium. 

There have been other attempts at analyzing microstrip antennas as open micro- 
strip cavity resonators for rectangular and circular patches (refs. 7 and 8). How- 
ever, these formulations do not include the patch excitation in the analytical model- 
ing, and the techniques cannot be extended to microstrip arrays. 

The purpose of the present work is to reformulate the problem to overcome the 
difficulties encountered in previous methods and to establish a foundation for ana- 
lyzing a variety of electromagnetic problems related to printed-circuit antennas. 

The problem is analyzed by first deriving a dyadic Green* s function which satisfies 
the boundary conditions for a unit current located in the plane of the microstrip 
patch. By weighting the Green's function with the electric current density and inte- 
grating over the patch, the radiated electromagnetic field is calculated at any point 
inside the dielectric. Coupled electric-field integral equations for the unknown 
patch current density are obtained by forcing the total tangential electric field on 
the patch to zero. With the proper basis and testing functions for the unknown cur- 
rent distribution, the integral equations are reduced to a matrix equation which may 
be solved for the patch current. The current distribution on the patch is then used 
to determine the properties of the microstrip antennas. 
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Some features of the present method are as follows: (1 ) it correctly accounts 

for the finite thickness of the dielectric substrate and the effect of surface waves? 
(2) it does not require high precision in the numerical computation of the matrix 
elements? (3) it allows improved convergence of the solution through proper choice of 
the basis and testing functions? (4) it can be extended to microstrip patches of any 
shape, provided one selects the proper basis functions? and (5) it is suitable for 
more involved problems of coupled microstrip antennas. 
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Eq( 9/ 4>) 

e .(e,4>) 
<t> 
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H i 
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magnetic vector potential 
dielectric substrate thickness 
electric-field vector 
amplitude of plane wave 

x component of electric field due to the incident plane wave (eq. (16)) 

y component of the electric field due to the incident plane wave 

0 component of E 

(j) component of E 

dyadic Green’s function 

bidimensional Fourier transform of G 

magnetic-field vector 

magnetic field due to the incident plane wave 

unit dyadic impulse current 

surface electric current density vector 
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j 

3x 

j y 

i x 
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k 
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rr 

bidimensional Fourier transform of J 
x component of j 
y component of j 

complex coefficient for x-directed current 
complex coefficient for y-directed current 
wave propagation constant in dielectric region 
wave propagation constant in free space 
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k x Fourier transform variable with respect to x 

ky Fourier transform variable with respect to y 

k z complex wave propagation constant in z direction 

M + 1 number of y-directed strips (see fig, 2) 
m,n m, n subdomain of field 

N + 1 number of x-directed strips (see fig, 2) 
p,q p,q subdomain of source 

r,0,4> coordinates of point in spherical coordinate system 

E-plane dimension of the rectangular patch 

(W /X ) resonant width of patch 

x' e o 

Wy H-plane dimension of the rectangular patch 

x, y,z rectangular Cartesian coordinates of field point 
x , ,y',z' rectangular Cartesian coordinates of source point 

AAA 

X, y,z unit vectors along x-, y-, and z-axis, respectively 
Z xx self -impedance of x-directed patch current (eq. (52) ) 


pq 

mn 

Z 

xx 


mutual impedance between x-directed currents on m,n and p,q subdomains 
(eq, (37)) 


pq 

mn 

Z 

xy 


mutual impedance between x- and y-directed currents on m, n and p,q 
subdomains (eq, (38) ) 


pq 

mn 

Z 

yx 


mutual impedance between y- and x-directed currents on m, n and p,q 
subdomains (eq, (39) ) 


pq 

mn 

Z 

yy 


mutual impedance between y-directed currents on m, n and p,q subdomains 
(eq, (40)) 


wave reflection coefficient 


Ax = W x /M + 1 

Ay = W y /N + 1 

A(W /X ) bandwidth for voltage standing-wave ratio (VSWR) less than 3 

x e 3 
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1 r°° 

ofx-x 1 ) delta function, — - I exp[j(x - x' )w] dw 

2u 


relative permittivity of dielectric substrate 
wavelength in dielectric substrate 
permeability of dielectric substrate 
to angular frequency, rad/sec 

Superscripts: 

I dielectric region 

II free space region 

I/II either region I or II 

THEORY 

General 

The geometry of the microstrip antenna is shown in figure 1. It consists of a 
thin conducting plate, or patch, in the form of a rectangle, a square, a circle, or 
other geometric shapes embedded in a dielectric substrate, with one surface being 
grounded by a perfectly conducting plane. Usually the microstrip antenna is excited 
either by a microstrip transmission line connected to one edge of the patch or by a 
coaxial probe with the center conductor extending through the ground plane and con- 
nected to the patch. However, in order to determine the resonant frequency, the 
bandwidth, and the radiation pattern, it can be assumed that the patch is excited by 
a plane wave polarized in the x-direction and with normal incidence on the patch, as 
shown in figure 1. This model gives results which are valid for the practical micro- 
strip antenna so long as the exciting mechanism does not appreciably alter the cur- 
rent distribution of the patch. For the dielectric substrate with a thickness less 
than 0.1 which is the case for most of the practical microstrip antennas, the 
dimensions of the exciting mechanism are much less than 1.0\ £ and, therefore, the 
feed connection produces only a small perturbation in the radiation characteristics 
of the antenna. 

The surface current induced on the patch controls the antenna properties such as 
radiation pattern, input impedance, resonant frequency, and bandwidth of the micro- 
strip antenna and, therefore, must be determined first. Let Jtx^y* ,z*) be the 
unknown surface electric current density present on the patch. For a very thin patch 
in the plane z = z', J has only x- and y-directed components. However, for the 
general case, it is assumed that J has all three components. The electromagnetic 
field due to J can be obtained from 

H(x,y,z) = — V x A(x,y,z) (1 ) 

r 


E(x, y, z ) = - ^{k 2 A(x,y,z) + V[V • A(x,y,z)]} (2) 

k 
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where the assumed time variation exp(jurt) has been suppressed. The magnetic vector 
potential A must satisfy the wave equation 


2-1 2-1 - 
V A (x,y,z) + k A (x,y,z) = -|iJ 

in the dielectric (region I) and must satisfy 


(3) 


V 2 A II (x / y / z) + k 2 A IJ (x,y,z) = 0 
o 


(4) 


in free space (region II). 

If G I ^ II (x / y / z f x' f y' f z' ) is the dyadic Green’s function for regions I or II 
for a unit impulse current source Kx'^'rZ 1 ) embedded in the dielectric substrate, 
then the solutions of equations (3) and (4) can be written in the form 


A 1 / 11 


(x,y,z,z' ) 


/ / J(x' ,y' ,z' ) 
Patch 


G I / II (x r y,z,x' , y* ,z' ) dx' 


dy 1 


(5) 


With the dyadic Green's function from the appendix substituted into equation (5) and 
the order of integration changed, the magnetic vector potential due to the patch 
current J(x',y f ) can be obtained in the form 


A 1 (x , y , z ) = 


/ /°° m j(k ,k ) • f 1 (k , k ,z) exp [ j (k x + k y) ] dk dk 

J —CO J —CO V \7 C V V \ 


=1 


(2 %y 


( 6 ) 


where 


j (k ,k ) = j j J(x*,y*) exp [-j (k x* + k y')] dx' dy ' 


x y 


Patch 


x y 


(7) 


and the components of the dyadic f are defined in the appendix. 
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The scattered field due to the patch current can be obtained from equation (2) 
in the following form: 


E I (x,y,z) = — ~ 3M - r I” [C (k ,k ,z) j (k ,k ) 
x ” .2 2 J - co L xx x y x x y 

\ ) JC 


+ C yx (k x ,ky, z) j y (k x ,k y )] exp(jk x x + jk y y) dk x dk y 


( 8 ) 


E X (x,y,z) = f P [C (k ,k ,z) 3 (k ,k ) 

y ' ,y ' ' o, 2 J L xy x y x x y 

y (2%) k 


+ C yy (k x'V Z> j y (k x' k y )] exp(jk x x + jk y y) ^x ^y 


(9) 




_I/ v -]U 

E (x,y,z) = 

Z / n \ Z 


(2 it) k *'- co 



k 2 + k 2 
x y 




f 1 (k ,k ,z) + f 1 (k ,k , z ) 
zx x y zy x y 


>] 


where 


5 f 1 (k ,k ,z) -.'j 

* 3 ” te y - [Vy'VV * k A°VV]J •*><»** + 3 V> ^y 

( 10 ) 


5f 


c „„, k „. V z, . (* 2 - ^)4 ♦ 


XX X 


( 11 ) 


5f 


C (k /k ,z) 

xy x y 


= -k k f + jk k 


zx 


A. i- I J 

x y xx x y oz 


(1 2 ) 


C (k ,k ,z) 
yx x y 


df 

-k k f + jk k ■ ~~r~~ 
x y yy x y oz 


(13) 
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( 14 ) 


C (k ,k ,z ) = (k 2 - k 2 )f I + jk 
yy x y y yy y 


5 f 1 

2 zy 


dz 


The total tangential electric field at the plane z = z' consists of the scat- 
tered field due to the patch current and the field due to the incident plane wave. 
The scattered field due to the patch current is given by equations (8) and (9), and 
the electric field at the patch location z = z' with the patch removed is obtained 
by using the laws of reflection and refraction to produce the following: 


l 

E = E 
x o| 


2j sin(kz') exp(jk Q d) 


^ eT cos(kd) + j sin(kd) 


(15) 


By setting the total tangential electric field over the patch equal to zero, the 
following integral equations are obtained: 


l 

E 

X 




(271) k 


2,2 C C, WV 


C (k ,k ,z' 
yx x y 


j (k ,k )] exp(jk x + jk y) dk dk 
yxy J x y x y 


(16) 


E = 


JO) 


(2ti) k 


o o / J [C (k , k , z 1 ) j (k , k ) 

2, 2 J J -<*> L xy x y x x y 


C (k , k , z 1 
yy x y 


j (k ,k )1 exp(jk x + jk y) dk dk 
yxy J x y x y 


(17) 


The integral equations (16) and (17) can be used to determine the patch current 
for a microstrip patch of any arbitrary shape. However, when one tries to solve 
equations (16) and (17) by the method of moments, it is difficult to arrive at a 
proper expansion function for the current distribution on the patch unless the geom- 
etry of the patch is an ordered shape such as a rectangle, a circle, or a square. It 
may be concluded that the application of equations (16) and (17) may be limited by 
complexity of the patch shape. 


Radiation Pattern 

The patch current obtained after solving equations (16) and (17) can be used to 
evaluate the radiation pattern. It is well known that the far field due to a radiat- 
ing aperture is directly related to the Fourier transform of the tangential electric 
field over the aperture (ref. 9) . Hence, the radiation pattern due to a microstrip 
patch embedded in a dielectric substrate with one surface grounded can be obtained 
from a knowledge of the Fourier transform of the electric field over the interface 
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between free space and the dielectric region, Ihe transformed electric field in the 
plane of z = d can be obtained from equations (8) and (9) as 


wv z) 


z ==d 


= 3-gijg r c (k ,k ) j (k ,k ) + c (k ,k ) j (k ,k 
|jxx x y J x x y yx x' y J y x ; 


■3 


e (k ,k f z ) 
y x y 


z=d 


-ju\)a 
, 2 


C 


C (k ,k ) j (k ,k ) + C (k ,k ) j (k ,k ) 

xy x y J x x y yy x y J y x y 


) 


where e x and e y are the Fourier transforms of and E^. The far field 

then given by (ref. 9) 


E. = [e cos 4> + e sin ^ exp(-jk r) 

0 ^ x y / 27tr o 


\ 


j* 


E f = fe cos 6 - e sin (J)] cos 0 w exp ( -jk 
4> l y x j 27ir * J 

With the following coordinate transformation 


r ) 

o 


k = k cos (J) sin 0 
x o 


k = k sin 4) sin 0 

y o 


the far-field radiation patterns become 


-k 0) exp(-jk r) 


E 0 (0,4>) = 


27irk 


-{t 


. . ,C (0, <t> ) cos <e + c (0 

2 M xx xy 


, <t> ) sin 4>j 


(6,4) ) 


C 


+ C (9,4)) cos 


yx 


<t> + C yy (9,4>) sin 4>J j y (9,4>)j 


E.(9,4>) = 

<P 


-k u exp(-jk r) 
o o 


2itrk 


I 


cos 0 / |c ( © / 4> ) c °s 4> - C (0f4>) sin (M j (0,4>) 
xy xx l x 


(1 8 ) 

(19) 
is 

( 20 ) 
( 21 ) 

(22a) 

(22b) 

( 23 ) 



r 


where the quantities C xx ( ^ # <t> ) r ^yx* 9 '^' and C yy (0,4>) are evaluated 

at z = d. Also, j x (0f<i>) and are Four ^- er transforms of the patch 

current with k and k_ r replaced by equations (22)* 
x y 

After considerable mathematical manipulation, the radiation patterns can be 
written in the form 


Eq( 9,<)>) = C o [j x (9,<}>) cos <j) + jy(9,4>) sin (j>](e r - sin^9)^|j 1^ 


sin ( D k - Z’) 


°k ' K 




d te {0 ' <},) 


Sin °k\ 2 

_ — - j (e - 1) sin 9 sin Z' 

D, / r k 


1 

1 

/sin D k \| 

d te < 9 '4>> 

d tm (0 '* ) 

{ )] 


(25) 


E^(9, <(>) = -C Q [j 9, 4>) sin <|) - j ^ ( 0 , c|) ) cos 4> ] ( cos 9) 


(f - ’) 


sin(D k - Z£) 


D k - K 


N TE (6 ' <t,) 


D te ( 9 , (J) ) 


^sin D, 


where 


(26) 


N (0,4)) = Vie - sin 2 0 cos(D - Z* ) + j cos 0 sin(D. - Z ' ) 
TE I r k k k k 


, $) = \[ e - sin 2 0 cos D, + j cos 0 sin D, 
TE T \ r k J fc 


D ( 0, ^) ) = e cos 9 cos D. + j \ e - sin 2 9 sin D, 
TM r k l|r ^ 


D. = k d \ e - sin 2 ( 
k o \ r 

Z ■ » k z 1 \| £ - sin^ 

k o V r 


( w ^o d | 

exp(-jk r) 
o 

\ 2 7l£ / 

k r 

\ r / 

o 
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Application to Rectangular Patch 


Consider a microstrip antenna in the form of a rectangular patch with respective 
E- and H-plane dimensions of W x and W , as shown in figure 2. The current distri- 
bution on the patch must be expanded in terms of eigenfunctions which satisfy the 
proper edge behavior such that the normal component of the patch current density is 
zero at the edge and the parallel component approaches infinity at the edge. In the 
classic eigenfunction solution of a problem involving edges, these edge conditions 
are explicitly used to select the eigenfunctions. However, it may be unnecessary 
here to use the exact eigenfunctions. Instead, one may expand the patch current 
density in terms of triangular, sinusoidal, or flat pulses to approximate the current 
distribution. Therefore, the rectangular patch is divided into rectangular subdo- 
mains which allow one to expand the current in terms of overlapping triangular pulses 
in the direction of the current and in terms of flat pulses orthogonal to the cur- 
rent. Hence, 


M N+1 


J (x 
x 


- y) -EE 

m=1 n=1 


i™ P U> Q ( y) 

x m n 


M+1 N 
V* :mn 

y 'm n 


V x ' y) ‘2^ 2^ i:r Q - (1 ° p - !y) 


m=1 n=1 


where 


P (x) = 1 + (x - x ) /Ax 
m m 

= 1 - (x - x ) /Ax 
m 

= 0 


( (x - Ax) < x < x ) 
m m 

fx < x < (x + Ax)) 
v m m ' 


Q^(y) = i (y < y < (y n - Ay)) 

= o 

Ax = W ^(M + 1 ) 

Ay = W / (N + 1 ) 

y 


(27a) 


(27b) 
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and iT an< ^ ^y 11 are t * le complex x- and y-directed current coefficients 

for the m, n subdomain. 

With equation (7)/ the Fourier transforms of the patch current from equa- 
tions (27) are 


where 


M N+1 


j (k ,k ) = V N j mn F (k ,k ) 
x x y /^i / ^ J x x r mn x y 


m=1 n=1 


M+1 N 


j (k ,k ) -V* j ™ 1 P (k ,k 

y x y / ^ / v J y y,mn x \ 


m=1 n=1 


(28a) 


(28b) 


J fY /*x +Ax 

11 j m P m (x ) {^(y) exp [-j (k^x + k y) ] dx dy 
v — Av •'x -Ay m x y 


F (k ,k ) 
x,mn x y 


y “Ay •'x -Ax 
n m 


= Ax Ay 


sin(k y Ay/2) 


k y Ay/ 2 


sin(k x Ax/2) 


k x Ax/2 


X exp[-jk x x m - jk y y n + jk (Ay/2)] 


(29) 


J r y +Ay /.x 

I ^n (x) p n (y) e3 ^p[*3< k x x + k „y)] dx d y 

v -Av •'x -Ax ^ 


F (k 
y ,mn x 


y -Ay •'x -Ax 
n m 


sin(k^ Ay/2) 

2 

sin(k x Ax/2) 

k y Ay/ 2 


k x Ax/2 


X exp[-jk X - jk y + jk (Ax/2)] 
xm yn x J 


(30) 


1 1 


Substituting equations (28) into equations (16) and (17) results in the following 
integral equations: 


E 1 « — 

x %2, 2 
(2n) k 


M N+1 

S E i“ r r c (k ,z') F (k ,k ) 

j v J —oo * — oo YV y V- . inn v \r 


— xx x y 
m =1 n =1 


x,mn x y 


M+1 N 


exp(jk x + ]k y) dk dk + Tl J™ P / C (k ,k , z ' ) 

x y x y " " J y ^ J -<» yx x v 

m =1 n =1 


yx x y 


x F (k ,k ) exp( jk x + jk y) dk dk 
y,mn xy ^ J x y x y 


(31) 


E"*" =- ^— 

y (2it) 2 k 2 


M N+1 _ 

Y Y j mn itoc /-oo 5 < k * k »z") F (k ,k ) 
^ ^ *'^x> ^xy x y x,mn x y 


m= 1 n= 1 


M+1 N ^ 

x exp(jk x + jk y) dk dk + Y' Y^ j mn /°° / C (k ,k , z' ) 

J x. y x y * J y ' — « oo * -oo yy x y 

m =1 n =1 


x F (k ,k ) exp( jk x + jk y) dk dk 
y,mn x y x y x 3 


(32) 


If the method of moments is used with testing functions identical to the expan- 
sion functions, equations (31) and (32) can be reduced to algebraic equations. Mul- 
tiplying equation (31) by P (x) Q (y) and integrating with respect to x and y 
gives P q 


M 

N+1 

pq 

M+1 

N 

pq 

= E 

E 

,mn mn 
j z 

X XX 

+ E 

E 

' mn mn 

j z 

y xy 

m =1 

n =1 


m =1 

n =1 


(33) 


Likewise, multiplying equation (3 2) by P (y) Q (x) 

q P 


and integrating gives 


M 

N+1 

pq M+1 

N 

pq 


E 

m =1 

E 

n =1 

, inn mn 
j Z + > 

J x yx 

m =1 

E 

n =1 

, mn mn 

j z 
y yy 

(34) 
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where 


E x,pq 



x +Ax 
P 

-Ax 

P 


E 1 P (x) 
x p 


Q (y) dx dy 
4 . 


E 

y/pq 



p q(y) 


Q p (x) dx dy 


pq 

mn io} ,«> 

Z - J J 5 (k ,k ,z') F (k ,k ) F (-k ,-k ) dk dk 

xx 2 2 -oo -oo xx x y x,mn x y x,pq x y x y 

(27i) k 


pq 

mn i a) r™ r°° 

Z = J J c (k , k , z ' ) F (k , k ) F ( -k , -k ) dk dk 

xy 2 2 - 00 -o° yx x y y,mn x y x,pq x y x y 

( 2tc) k 

pq 

mn i co r °° 

Z - — J J C (k , k , z 1 ) F (k , k ) F (-k ,-k ) dk dk 
yx , 2 2 xy x y x,mn x y y,pq x y x y 

( 2 tt) k 


pq 


mn 


Z 

yy 


, 2 2 
( 2 ti) k 



/ C (k ,k ) F (k ,k ) F (-k ,-k ) dk dk 
yy x y y,mn x y y,pq x y x y 


By varying p = 1, 2, M and q = 1, 2, . .., N + 1 in equation (33) and by 

varying p = 1, 2, . .., M + 1 and q = 1, 2, ..., N in equation (34), 

M(N +1) + (M + 1)N equations with M(N +1) + (M + 1)N unknowns are obtained. 
Hence, equations (33) and (34) can be written in the following matrix form, where 
dimensions of each submatrix are given below the submatrix: 


- 

— 


— 

“ pq* 

1 

1 

~ pq' 




r~ - 



r i 



mn 


mn 




.-mn 



E 1 



Z 

1 

Z 




i 



_ X,pql 



XX 

1 

- 




L x J 


M(N+1 ) xl 


M(N+1)xm(N+1) J M(N+1)x(m+1)N 


M(N+1 ) x 1 





" pq~ 


-pq" 







r 1 



mn 

1 

mn 




’-mn 



E 1 



Z 

1 

Z 




i 



Ly'Mj 



_ v*. 

1 

1 

yy 




_ y 


< 

M+1 )NX1 


(M+1 )N X M( 

N+1) I (M+1) 

NX (M+1 )N 


(M+1 )NX1 

_ 







-J 


L. 


— j 


(35) 


(36) 


(37) 


(38) 


(39) 


(40) 


the 


(41) 
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k * k 8 sin a, the gen- 
y 


With the variables changed such that k x = k^p cos a and 
eralized impedance elements in equation (41) are given by 


pq 

mn 

Z = C 
xx 


r r 


b (a, p) 


sin| — k 8 Ax cos a] 
\2 o )_ 


-k 8 Ax cos a 
2 o 


p=0 a=0 

x [(e - p 2 cos 2 a) 5 (p) - j (e - 1)p 2 cos 2 a C_(P)] 

r I r z 

X cos [k p(x - X ) cos a] cos[k p(y - y ) sin a]p dp da 
L o^ m p J L o r J n q 


pq 

mn 

Z = C 

xy 


J ~n/2 
n a a 


b(a,p) 


P=0 a=0 


— 

/i \ 

sinl 

f 1 \ 

— k p Ax cos a 


\,2 o J 

1 


— 

k p Ax cos a 

_ 2 

o 


sinfl k p Ay sin a) 
\2 o / 


— k P Ay sin a 
2 o 


x [5,0) + j(e - D S 2 (P)](3 2 sin a cos a 

in [v( X m - x p + f 1 ) COS “] Sin [ k o P ( y n " y q + f) Sin a ] P dp 


x sm 


da 


pq tnn 

mn pq 

Z = Z 
yx xy 


pq 

mn 

Z = C 

yy 


J r m [K/2 

n r\ ~ r\ 


b(a f P) 


sinf— k P Ay sin a) 

\2 o ) 


1 

— k p Ay sin a 
2 o 


'p=0 a=0 

X [(e r - p 2 Sin 2 a) C, <P> - j (e r - 1 )p 2 sin 2 a ? 2 (p)] 

x sinfk P(x - x ) cos a] sin[k P(y - y ) sin alp dp da 
L o K m p J L o r n q 


where 


C = 3 120 (k z’ ) (k Ax) (k Ay) 
ue o o o 

r 


1 4 


(42) 


(43) 


(44) 


(45) 


( 46 ) 




and the integration over the k x ,k -plane has been reduced to an integration over one 
quadrant by using the even and oda properties of the integrand. 


The impedance elements are computed by using numerical integration techniques 
and properly accounting for the contribution of surface wave poles. The unknown 
complex current coefficients and are then determined from equation (41) 

by using matrix inversion. 

Patch Current Distribution 

The patch current distribution can be obtained by solving the matrix equa- 
tion (41). The matrix elements are numerically computed by using the Gauss 
quadrature numerical integration technique. For accurate computation, the upper 
limit on the (3 integration must be terminated ait a sufficiently large value. It 
was found numerically that any increase in the upper limit above 600 does not change 
the values of the matrix elements appreciably. 

Without loss of generality, the left-hand side of equation (41) may be assumed 
to be unity. Therefore, the numerical results for patch current distribution given 
in this report are normalized to the total incident electric field on the patch, that 
is, E x given in equation (15). 
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Although the problem of the rectangular microstrip antenna is formulated cor- 
rectly in terms of coupled integral equations, the approximation to the patch current 
distribution is introduced when one tries to solve the integral equations by the 
method of moments or by Galerkin's method. The numerical solution of equation (41) 
is expected to converge when the number of subdomains approaches infinity. In prac- 
tice, however, because of computer storage and time limitations, it is desirable to 
obtain an answer with as few subdomains as possible. It is, therefore, important to 
know the proper number of subdomains to obtain accurate values of the patch current 
distribution. 

For a 1.0\ square patch (W v = W ) in free space and illuminated by a plane 

G ^ Jt 

wave at normal incidence, the patch current converges to the correct value when 
M > 7 and also N > 1 (ref. 10). The present problem becomes identical to the one 
in reference 10 when d = 00 and = 1 . It is expected that the patch current 

obtained after solving equation (41) for M > 7 and N > 7 will oscillate about the 
converged value of the patch current in reference 10 as the substrate thickness d 
is varied. The patch current at the center of the patch obtained after solving equa- 
tion (41) for M = N = 7 is plotted in figure 3 as a function of the substrate 
thickness for a dielectric constant of unity. As expected, both the real and imagi- 
nary parts of the patch current indeed oscillate about the free space value. For 
further calculations, the number of subdomains is selected such that M = N = 11. 

For further verification of the present method, it is important to check the 
results obtained with the present method with earlier known results. The transverse- 
electric-strip (TE-strip) problem is a limiting case of the rectangular patch under 
the condition Wy ■* °°. The patch current at the center, therefore, must asymptoti- 
cally approach the TE-strip current (refs. 11 and 12) as Wy approaches infinity. 
Figure 4 shows the results are as ejected. 

After checking the validity of the present formulation, equation (41) is solved 
for a square patch. For x-polarized plane wave excitation, the magnitude of the 
y-directed current is very small compared with the magnitude of the x-directed 
current (typically 25 to 30 dB less); therefore, Jy may be neglected. It is 
advantageous to study the current distribution on the patch, J x (x,y), which may 
allow one to express the current in terms of more appropriate basis functions. The 
current distribution J x (x,y), obtained after solving equation (41), is plotted in 
figure 5 as a function of x at y = Wy/2 and as a f unct i° n of y at x = W^/2. 

It is clear from figure 5 that J x (x,y; is almost sinusoidal in the direction of the 
current and is uniform, except at the edges, along the direction orthogonal to the 
current; therefore, it is more appropriate and advantageous to express the patch 
current density in the following form: 


J (x,y) = i cos(ttx/W ) (50) 

x J x x 


The assumption of uniform distribution in the y direction gives a value of J 
slightly higher than the true value. With the above current distribution, the number 
of unknowns to be evaluated is reduced from M(N + 1) + (M + 1)N to 1 . The matrix 
equation (41) then reduces to 


x 


= Z j 

XX X 


(51) 
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where 


r°° rV 2 

sinf— k BW sin a) 

2 

2n cos(“ k Bw cos a) 

1 

\2 o y / 


\2 o x / 

(3=0 *a=0 

1 k Bw sin a 

2 cT y 


2 2 
ix - (k BW cos a) 

o r X 


[<e - P 2 


2 x 

cos a) 


5 1 <P> - 


j(e - 1)p' 


2 

cos a 


5 2 <p>]p dp da 


( 52 ) 


and C, ^(p), and £ 2 (f3) are given equations (46), (48), and (49). 


RESULTS 

Patch Resonance 

It is informative to plot the real and imaginary parts of / as a function of 
the E-plane dimension, as shown in figure 6. The plot has the characteristic reso- 
nant behavior of a tuned circuit and can be used to predict the resonant width (i.e., 
the resonant frequency) of a rectangular microstrip antenna. The resonant frequency 
can be defined as the frequency at which the imaginary part of J x is zero. With 
the above criteria, the resonant frequency of a rectangular patch is determined and 
plotted as a function of the H-plane dimension in figure 7 along with available theo- 
retical and experimental results from reference 4. It should be noted that the theo- 
retical prediction of resonant frequency using the cavity model (ref. 4) requires a 
semiempirical correction. However, the present method of analysis, without any cor- 
rection, predicts the resonant frequency within the tolerance limits of the dielec- 
tric constant (e = 2.62 ± 0.02). 
r 


Resonant Width 

With the resonant criteria defined in the previous section, the resonant width 
of a rectangular microstrip patch is computed as a function of dielectric substrate 
thickness for the aspect ratios W /W^ of 0.5, 1.0, and 1.5 and for dielectric con- 
stants e of 1.0, 2.5, 6.0, and T0.0. The values are presented in figures 8 
to 11. Tfiese calculations show that the resonant width (normalized to the dielectric 
wavelength) of a rectangular microstrip patch antenna generally increases with 
increases in the dielectric constant of the substrate and decreases with increases in 
the substrate thickness. 

The resonant width of a rectangular patch is also plotted in figure 1 2 as a 
function of the aspect ratio W y/ W x * An increase in the aspect ratio causes a 
decrease in the resonant width of the patch. Under the limiting case of W y/W x 
being large, the resonant width of the rectangular patch will approach asymptotically 
the resonant width of the TE-strip (ref. 12). 

It is informative to plot the resonant width of the rectangular patch for a 
thick dielectric substrate, as shown in figure 13. These calculations show the reso- 
nant width decreases approximately linearly with increases in substrate thickness 
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for d < 0.06X • For d > 0.06X , the curve is no longer linear and scaling can no 
longer be used to estimate the resonant frequency of a microstrip antenna# 

In some applications it is desirable to cover the microstrip patch with a 
dielectric coating to protect the antenna from environmental conditions# The reso- 
nant width of a rectangular microstrip patch for various dielectric coating thick- 
nesses is computed and is presented in figure 14. Hiese calculations indicate that 
the resonant width decreases with an increase in the thickness of the coating# 


Bandwidth 

Hie bandwidth of microstrip antennas may be defined as the band of frequencies 
over which the input voltage standing-wave ratio (VSWR) in the feed line is less than 
a specified value {usually 2 or 3). Since the input impedance of a microstrip 
antenna is proportional to the patch current, the input impedance variation with 
frequency will be the same as that of the patch current. Hence, the input reflection 
coefficient variation with frequency can be calculated from a knowledge of the patch 
current density# Under the assumption of a unit voltage source and a perfect imped- 
ance match at resonance, the reflection coefficient T can be calculated from 


~ -^x,r 
ix + ^x,r 


(53) 


where j is the amplitude of the patch current density at resonance. Hie magni- 
tude of t&e input reflection coefficient thus calculated is shown in figure 1 5 as a 
function of the patch size for a square patch. For the input VSWR in the feed line 
to be less than 3, the percent bandwidth may be obtained by 


Percent bandwidth 


100 A(w x A £ ) 3 

(W /\ ) 

X £ O 


(54) 


where A(W x /\ £ ) 3 an< * (W x A e ) 0 are figure 15. The percent bandwidth for 

a VSWR less than 2 can be determined in a like manner. Hie percent bandwidth thus 
calculated is compared with measured data in figure 16 for probe-fed square patches. 
The apparent increase in the measured percent bandwidth may be due to feed-connection 
insertion loss and copper loss in the experimental models. Hie calculated percent 
bandwidths for rectangular patches are presented in figures 17 to 22 for aspect 
ratios W y/W x of 0.5, 1.0, and 1.5 and dielectric constants e ^ of 1.0, 2.5, and 
6.0. These calculations show an increase in bandwidth with increases in dielectric 
substrate thickness, which is consistent with previous measured results (ref. 13). 

The percent bandwidth of a square microstrip patch with a thicker dielectric 
substrate is presented in figure 23. This calculation indicates that a bandwidth up 
to 20 percent may be possible by constructing the antenna on a thick dielectric. 

Figure 24 shows the effect of a dielectric coating on the percent bandwidth of a 
square patch. Hie calculations show no appreciable change in percent bandwidth when 
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the cover thickness is less than 0.04X • However, for thicker dielectric coatings 
(greater than 0.04\ £ ), there is an increase in the bandwidth of about 1 percent. 

Radiation Pattern 

With the x-directed current distribution given in equation (51), the radiation 
patterns of a rectangular micros trip antenna are computed from equations (25) 
and (26). The results are presented in figures 25 and 26 along with experimental 
data (ref. 4). There is a close agreement between the calculated and measured 
radiation patterns. 


CONCLUSION 

The dyadic Green's function for a current source located in a dielectric with 
one surface grounded has been derived. The dyadic Green's function was used to 
arrive at a set of coupled electric-field integral equations for a microstrip patch 
of arbitrary shape. As an example of the applicability to microstrip antennas, the 
integral equations were solved for a rectangular patch using Galerkin's method. It 
is shewn that the patch current can be used to accurately predict the resonant fre- 
quency, the bandwidth, and the radiation patterns of microstrip antennas. Data for a 
range of design parameters are presented for the rectangular microstrip antenna. 


Langley Research Center 

National Aeronautics and Space Administration 
Hampton, VA 23665 
January 25, 1984 
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APPENDIX 


DERIVATION OF DYADIC GREEN'S FUNCTION 

The expression in equation (5) must satisfy the wave equations (3) and (4); 
therefore. 


V 2 G I (x,y,z,x' ,y' ,z' ) + k 2 G 1 (x, y, z,x' , y' , z ’ ) = -|il d(x-x’) 6(y-y') S(z-z') 


(A1 ) 


V 2 G Z1 (x,y,z,x’ ,y' ,z' ) + k 2 G II (x,y,z,x',y',z' ) = 0 

o 


(A2) 


It is clear from equations (A1 ) and (A2) that the Green's function is the electric 
vector potential due to a unit impulse current located at the point (x^y^z 1 )* 

It is advantageous and less tedious to obtain the solution of equations (A1 ) 
and ( A2) in the spectral domain* Therefore, taking the bidimensional Fourier 
transform of both sides of equations (A1 ) and (A2) and using Sommerf eld's radiation 
condition (ref. 14) yields 


aV(k x ,k ,z,x',y',z') ml 

+ (k ) g (k ,k , z,x 1 , y' , z 1 ) 

bz y 


= -|i.I 6(z-z') exp[-j (k^x' +k^y')] 


(A3) 


& 9 (k ,k ,z,x' ,y' ,z' ) 

+ (k ) g (k ,k , z,x ' , y' , z ' ) = 0 

5z 2 x Y 


(A4) 


where 


g I//II (k x ,k y ,z,x' ,y' ,z' ) = j_ m G 1 ^ I1 (x, y, z,x ' , y' , z ' ) ejq>[-j(k v X + k^y) ] dx dy 


x y 




(A5) 

k 1 = |k 2 e - (k 2 + k 2 ) 
)| o r ^ x y 

(k 2 e > (k 2 + k 2 ) ) 
v o r v x y 

( A6a) 

= -j i (k 2 + k 2 ) - k 2 e 
('■ x y ' or 

(k 2 e < (k 2 + k 2 ) ) 
v o r v x y j ' 

(A6b) 
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k 11 = ,|k 2 - (k 2 + k 2 ) 
I o v x Y 


(k 2 > (k 2 + k 2 )) 
v o v x y* 1 


(A7a) 


• - J \I( k x + k y) - k , 


(k 2 < (k 2 + k 2 )) 

k o v x y J ' 


(A7b) 


With a shortened notation g 1 ^ 11 , g I ^ II (k ,k ,z,x' ,y' ,z* ) may be expressed as 

x y 


f I/II , gl/llj + 5 X/H- + 5 I/HS 


(S 8 ) 


where 


9 1 / 11 - + g ^ 11 ? + 9V 11 * 


3 XX 


(A9) 


?y iz - 9ic zi i * + g yy-z 


xy 


"yy 


3 zy 


(A1 0) 


-I/I I = I/I I- + I/Ii; + I/II^ 
y z y xz ^yz y zz 


(A1 1) 


The scalar components of g 1 ^ 11 satisfy the following wave equations: 


jv 2 I 

d g 

xx 

^ 2 
5z 


+ (k 1 ) 2 g^ x = -p 6 (z-z') exp[-j(k x x’ + k^y')] 


(A1 2) 


*2 1 

!!zZ . * 


9 z 


+ (k ) g yy = -p 6 (z-z ’ ) exp[-j(k x x’ + k^y * ) ] 


(A1 3) 


a 2 1 

6 g zz fl I >2 I 


5 z 


+ (k ) g zz = -p 6 (z-z • ) exp[-j(k x x’ + k^')] 


(A14) 
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Other components satisfy the following equations for regions I and II: 


5z 2 


oo 


f 1 = 0 


b 2 f IX 

az 2 


+ (k 11 ) 2 f 11 = o 


(A1 5) 


(A16) 


where f x denotes the right-hand scalar components of equations (A9) to (All) 


(excluding g 


xx' 


g yy' 


and 


^zz 


for region I, and f 


II 


denotes the right-hand 


scalar conponents of equations (A9) to (All) for region II« 


The solution of equation (A12) has two parts, complementary and particular, 
complementary solution can be written in the form 


The 


I 

g 

xxc 


C + exp(-jk I z) + C exp(jk I z) 

XX XX 


(A17) 


where C xx and C” x are unknown constants. For the particular solution, the 
Fourier transform of both sides of equation (A1 2) is taken with respect to z and 
yields 


g (k ,k ,k ,x 1 , y* , z ' ) 
xxp x y z 


H exp[-j(k x x’ + k y y* + k z z 1 ) ] 

** o -- 1 ; 2 


(A18) 


Taking the inverse transform with respect to z yields 


g xxp (k x'V Z ' X, ' Y ’' Z,) = 


p exp[-j (k^x* 


+ k y’ )] 
y 


271 


r 


ex p[jk z (z - 


)] 


k z 2 - (*V 


dk 


(A19) 


The integrand in equation (A19) has poles at k z = ±k T ; therefore, equation (A19) can 
be evaluated through contour integration in the complex domain. Hence, 


g* = exp[-j (k X* + k y ' ) ] exp[-jk I |(z 

xxp 2k r x y 


Z ’)|] 


(A20) 
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The complete solution of equation (A1 2) can be obtained from equations (A1 7) 
and ( A20) as 


- exp[-j (k x' + k y ' )] exp[-jk I |(z - z' )| ] 

2k x T 


(A21) 


Likewise, the solutions of equations (A13), (A1 4) , (A15), and (A1 6) are obtained in 
the following forms: 


g 1 = C + exp( -jk X z ) + C exp(jk I z) 

yy yy yy 


- exp [-j (k v x' + k^y ' ) ] ejq>L-j kI U z ~ z ' ) I ] 


2k 


(A22) 


g zz = C zz ex P(-3 kI z) + c" z exp( jk*z ) 


- e^p[-j (k x* + k v' ) ] exp[ -jk 1 1 (z - z • ) | ] 

2k Y 


(A23) 


f 1 = 


C + exp(-jk I z) + C 


e3q>(jk z) 


(A24) 


,,11 II x 

f = D exp ( -jk z) 


( A25) 


where C + and C~ are the coefficients of g 1 excluding those given in equa- 
tions (A21) to (A23), and D + are the coefficients of g* 1 . 

The unknown coef f icients^i ^equations (A21) to (A25) are evaluated by applying 
the boundary conditions on g ' • These boundary conditions are obtained by equat- 

ing to zero the tangential electric field at the plane z = 0 and enforcing the 
continuity of tangential electric and magnetic fields at the plane z = d. The 
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Fourier transforms of the electromagnetic fields in regions 3 
through the following relations: 


h I/n (k ,k ,z,x' ,y' ,z' ) - 1 (k x g 1 ^ 11 ) 
x y 


e Iy/II (k ,k ,z,x' ,y’ ,z' ) 
x y 


■^(K x K 


x 5VII; 


where 


K = x(jk x ) + y(jk y ) + z 


(fe) 


Subjecting the field obtained from equations (A26) and (A27) 
tions gives, for zero tangential electric field at z = 0, 


I I I I 1 1 

g = g = g = g = g =g =0 
~xx xy ^ xz ^ yx y y z 


dg 1 dg 1 Qg 1 
zx zy _ ^zz _ 


= 0 


dz dz dz 

Continuity of the tangential magnetic field at z = d yields 


A 1 7K 11 

dg dg 

xx xx 


dz 

dz 

„ 1 

„ II 

5g 

dg 

xy 

xy 

dz 

dz 

^ I 

^ II 

dg 

xz 

dg 

~xz 

dz 

dz 

^ I 

^ II 

dg 


yx 

yx 


dz dz 


and II may be obtained 


(A26) 


(A27) 


(A28) 


to the boundary condi - 


( A29) 


(A30) 


(A3 1 ) 


(A3 2) 


(A3 3) 


(A34) 
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dg‘ 


yy _ 


5g‘ 


ii 

yy 


bz 


bz 


e>g 


yz 




II 

yz 


5z 


bz 


bg; 


zx 


&g; 


ii 

zx 


bz 


bz 


Sg; 


zy 


^g; 


II 

zy 


bz 


bz 


&g‘ 


zz 


^g; 


II 

zz 


bz 


bz 


I II 

g = g 

zx ^zx 


I ii 

g = g 
zy zy 


I II 

g = g 
^zz zz 


Continuity of the tangential electric field at z = d yields 


(A3 5) 
(A3 6) 
(A3 7) 
(A38) 
(A3 9) 
(A40) 
(A41 ) 
(A42) 


I II 

g = g 

XX XX 


I II 

g = g 

xy xy 


I II 

g = g 

xz ^xz 


i ii 

g = g 

yx yx 


g = g 

yy yy 


(A43) 


(A44) 


(A45) 


(A46) 


(A47) 
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I II 

g = g 

yz yz 


k x (e r " 1)g ^y +k y (e r ' 1)g yy = j l e i 


k (e - Dg 1 + k (e - l)g"„ = jle 
y r ^yx x r xx \ r 


A 11 

dg 

zy 

ag 1 \ 

zy 

dz 

dz / 

a 11 

dg 

zx 

dg 1 \ 

zx ] 

dz 

dz / 

* II 
dg 
zz 

dg 1 \ 

zz 1 

dz 

dz / 


Subjecting equations (A21) to (A25) to the conditions given by equations (A29) 
to ( A51 ) and solving the set of simultaneous equations for C and C results in 


g xx = f xx exp[ - j(k x X ' + V' )] 


g yy = f yy exp [- j(k x x ' + V'^ 
g zz = f zz e ^[-i (k x x ’ + V'^ 


g zx = f zx k x exp[ " j(k x x ’ + V' )] 


I 


zy zy y 

where, for (z - z') < 0, 


= f Z y k y + kyY 1 )] 


I _ I _ k 1 cos k 1 (d - z * ) + jk 11 sin k 1 (d - z* ) sin(k I z) 
xx ^ k 1 cosfk^d) + jk 11 sin(k I d) _ _ k*z 


[k 1 cos k I (d - z* ) + jk II e sin k 1 (d - z* ) 

3 M' . _ . r 


zz . I 
k 


k I:t e cos(k I d) + jk 1 sin(k I d) 
r 


cos(k z) 
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or, for (z - z 1 ) >0, 


f 1 

xx yy 


= -V*\ 


sin k (z - z 1 ) 


k x z 


Mz 


k 1 cos k 1 (d - z') + jk 11 sin k 1 (d 


- z 1 ) 


k 1 cos(k I d) + jk 11 sinCk^) 


sin(k z) 


k*z 


f = -uz 

17 17 1 


zz 


sin k (z - z* ) 


k X z 


JJi 

k 1 


k 1 cos k^d - z') + jk l:t e r sin k 1 (d - z* ) 


k II e^ cos(k I d) + jk 1 sintk'^d) 


{cos(k I z)J 


Also, 


^(e r - 1) sin(k I z l ) cos(k I z) 


f 1 = 

ZX Zy [k 1 cos(k I d) + jk 11 sin(k I d) ] [k I]C e cos(k T d) + jk 1 sinCk^)] 


The other components of g 1 are zero. Likewise, the components of g 11 are 
by 


g xx = f xx e *p[-3< k x *' + kyY’)] exp[ -jk 11 (z - d) ] 


„ II, 


g yy = f yy e *p["3 (k x *' + k yY'>] exp[-jk lx (z - d)] 


„ II, 


g zz “ f zz e ^p[-3( k x x ' + V’)] exp[-jk X1 (z - d)] 


„ II, 


g zx = f zx k x e3{ P[-i (k x x ' + V’)] exp[-jk IZ (z - d)] 


g ly = f zy k y e3 5>[~3 < k x x ’ + kyY')] exp[-jk i:E (z - d) ] 


(A59) 


(A60) 


(A61) 


given 


(A62) 


(A63) 


(A64) 


(A65) 


(A66) 
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where 


f 11 = f 11 
xx yy 


[j. sin(k I z ' ) 


k 1 cos(k I d) + jk 11 sin(k I d) 


(A67) 


f 11 = -\i(d - z' ) 
zz 


sin k (d - z 1 ) 


k (d - z 1 ) 


, I 


k cos 


k 1 (d - z ’ ) + jk II £^ sin k 1 (d - z 1 ) 


k 11 ^ cos(k I d) + jk 1 sin(k I d) 


Jcostk^jJ 


(A68) 


f 11 = f 11 


|j,(e - 1) sin(k I z l ) cos(k X d) 

r 


Zy [k 1 cos(k I d) + jk 11 sin(k I d) ] [k i:t £ cos(k I d) + jk 1 sin(k I d) ] 


(A69) 


Ihe other components of g 11 are zero. 

With the inverse Fourier transform, the dyadic Green’s function for regions I 
and II is obtained from 


/ \ 2 - 

G 1 (x, y,z,x ' , y' , z' ) = g 1 ^ ,k y , z,x ' , y* , z • ) exp[j(k x x + k y y)] dk x dk^ 


(A70) 


G II (x,y,z,x\y',z’ ) = (^j / P a g II (k x ,k y , z,x ’ , y’ , z ’ ) exp[j(k x x + k y y) ] 


. dk dk 

y J x y 


(A71) 
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Plane wave 



Region II 
(z > d) 



Figure 1.- Geometry of microstrip antenna. 
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Figure 4.- Comparison of rectangular patch current with TE-strip current 
W = 0.5\ ; d = z ' = 0.02X ; e = 2.5. 









Figure 7.- Resonant frequency of rectangular patch versus H-plane dimension 
W = 5.95 cm; d = z’ = 0.159 cm; e = 2.62 + 0.02 • 
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Figure 1 8 .— Calculated percent bandwidth for rectangular patch for 

VSWR <3; z' = d. 
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Figure 20.- Calculated percent bandwidth for square patch for 

VSWR <3. z' = d. 
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Figure 23.- Calculated percent bandwidth for square patch with thick dielectric 
substrate for VSWR <3. z* = d? e = 2.5. 
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Figure 24.- Calculated percent bandwidth for dielectric-coated square patch for 

VSWR <3. z’ = 0.0175X ; e = 2.5. 
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• Calculated 


(a) H-plane pattern. <j> = it/2. 

Figure 25.- Radiation pattern for 1187-MHz rectangular-patch antenna. 

W = 7.6 cm; W = 11.43 cm; z' = d = 0.158 cm; £ = 2.62. 
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— Measured (ref. 4) 
• Calculated 


(b) E-plane pattern. $ = 0. 


Figure 25.- Concluded 
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• Calculated 


(a) H-plane pattern. = t/2. 

Figure 26.- Radiation pattern for 804-MHz rectangular-patch antenna. 

W = 11.43 cm; W = 7.6 cm; z' = d = 0.158 cm; e = 2.62. 
x y r 
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— Measured (ref. 4) 
• Calculated 


(b) E-plane pattern. 4> = 
Figure 26.- Concluded. 
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